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Abstract 

We calculate the n + K + scattering length in fully-dynamical lattice QCD with domain- wall valence 
quarks on MILC lattices with rooted staggered sea-quarks at a lattice spacing of b = 0.125 fm, 
lattice spatial size of L = 2.5 fm and at pion masses of m^ ~ 290,350,490 and 600 MeV. The 
lattice data, analyzed at next-to-leading order in chiral perturbation theory, allows an extraction of 
the full irK scattering amplitude at threshold. Extrapolating to the physical point gives m^ a 3 / 2 = 
-0.0574±0.0016l&oo58 and m n a 1/2 = 0.1725±0.0017±g;g?§| for the I = 3/2 and I = 1/2 scattering 
lengths, respectively, where the first error is statistical and the second error is an estimate of the 
systematic due to truncation of the chiral expansion. 
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I. INTRODUCTION 



In hadronic atoms, nature has provided a relatively clean environment in which to explore 
the low-energy interactions of charged hadrons. The electromagnetic interaction allows for 
oppositely-charged, long-lived hadrons to form Coulomb bound states. The locations of the 
energy-levels of these systems are perturbed by the strong interactions, while the lifetimes 
of the ground states are dictated by the strong interactions that couple the charged hadrons 
to lighter neutral ones. 

Theoretically, the simplest hadronic atom to understand is composed of two pions: 7r + 7r~. 
Precision experiments have been performed — and are ongoing — to measure the lifetimes and 
energy-levels of such atoms In the isospin limit, Bose statistics dictates that two pions 
interacting in an s-wave can be in either an isospin-0 or isospin-2 state. By measuring the 
decay width and energy levels of pionium, the 1 = and 1 = 2 strong-interaction scattering 
lengths can be isolated. While the difference between energy levels, and hence deviations 
from the Coulomb spectrum, are relatively straightforward to measure, it is somewhat more 
challenging to determine the lifetime of these atoms. Recently, the DIRAC collaboration 0] 
at CERN has measured the lifetime to be r^+^- = 2.91^o;g 2 x 10~ 15 s, with the dominant 
decay mode being tt + tt~ — > 2tt°. On the theoretical side, progress in lattice QCD has been 
quite rapid, with a recent fully-dynamical calculation of the / = 2 tttt scattering length 
at pion masses between m n ~ 290 MeV and 500 MeV |2j. When combined with two-flavor 
chiral perturbation theory, a prediction of the scattering length at the physical point is found 
to have an uncertainty that is somewhat smaller than that from experiment. An up-to-date 
discussion of the status of ^-interactions can be found in Ref. 0]. 

Studying the low-energy interactions between kaons and pions with tt~K + bound-states 
allows for an explicit exploration of the three-flavor structure of low-energy hadronic in- 
teractions, an aspect that is not directly probed in irn scattering. Experiments have been 
proposed by the DIRAC collaboration 4] to study tiK atoms at CERN, J-PARC and GSI, 
the results of which would provide direct measurements or constraints on combinations of 
the scattering lengths. In the isospin limit, there are two isospin channels available to the 
tiK system, I = \ and I — %■ The width of a ir~K + atom depends upon the difference be- 
tween scattering lengths in the two channels, T ~ (ai/ 2 — CL3/2) 2 , (where ai/ 2 and a 3 / 2 are the 
I = \ and I = I scattering lengths, respectively) while the shift of the ground-state depends 
upon a different combination, AE ~ 2ai/ 2 + a 3 / 2 . Recently, the Roy-Steiner equations 
(analyticity, unitarity and crossing-symmetry) have been used to extrapolate high-energy 
ttK data down to threshold 0], where it is found that 

m n (ai/2 - a 3 / 2 ) = 0.269 ± 0.015 , m n (a 1/2 + 2a 3/2 ) = 0.134 ±0.037 , (1) 

which can be decomposed to m n ai/ 2 = 0.224 ± 0.022 and m 7T a 3 /2 — —0.0448 ± 0.0077. 
In addition, three-flavor chiral perturbation theory (xPT) has been used to predict these 
scattering lengths out to next-to-next-to- leading order (NNLO) in the chiral expansion. At 
NLO 0,0,3, 

m n (ai/2 - a 3/2 ) = 0.238 ± 0.002 , m n (a 1/2 + 2a 3/2 ) = 0.097 ± 0.047 , (2) 

while at NNLO 9] ^7^1/2 — 0.220 and = —0.047 1 . One must be cautious in 

assessing the uncertainties in these theoretical calculations, as one can only make estimates 



1 At tree-level, Weinberg |l0j determined that m^a^ = 0.137 and = —0.0687. 
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based on power-counting for the contribution of higher-order terms in the chiral expansion. 
There has been one determination of the tt + K + scattering length in quenched QCD 
however the chiral extrapolation was restricted to tree level. 

It is worth mentioning a novel motivation for accurate determinations of meson-meson 
scattering from lattice QCD calculations. Recent work has identified in a model-independent 
way the lowest-lying resonance in QCD which appears in nil scattering [12]. Crucial to this 
development has been the accurate determination of the low-energy tttt scattering amplitude, 
including the recent lattice QCD determination of the 1 = 2 scattering length [2j. A similar 
analysis has very recently been carried out for nK scattering in the I = | s-wave in order 
to determine the lowest-lying strange resonance [l3T ]. Improved accuracy in the low-energy 
tcK scattering amplitude should be welcome to this endeavor. 

In this work we present the results of a fully-dynamical lattice QCD calculation of n + K + 
scattering. By calculating the and m# dependence of the ir + K + (I = |) scattering 
length, we are able to provide a determination of both the / = | and I = \ scattering lengths 
at the physical point. We have performed a hybrid mixed- action calculation with domain- 
wall valence quarks tuned to the staggered sea-quark masses of the MILC configurations. As 
the computer resources do not presently exist to perform such calculations at or very near 
the physical value of the light-quark masses, these are performed at pion masses between 
m-n ~ 290 MeV and ~ 600 MeV. These results are combined with calculations in continuum 
three-flavor %PT to extrapolate to the physical point. 



II. FINITE- VOLUME CALCULATION OF SCATTERING AMPLITUDES 

The s-wave scattering amplitude for two particles below inelastic thresholds can be deter- 
mined using Luscher's method 3|> which entails a measurement of one or more energy 



levels of the two-particle system in a finite volume. For two particles with masses m,\ and 
m2 in an s-wave, with zero total three momentum, and in a finite volume, the difference 
between the energy levels and those of two non-interacting particles can be related to the 
inverse scattering amplitude via the eigenvalue equation 
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poot*(p) = 4 S (^) ■ ( 3 ) 
where 5(p) is the elastic-scattering phase shift, and the regulated three-dimensional sum is 

|j|<A i 

S M = E |Tp— 3 - 47rA • ( 4 ) 
j IJI / 

The sum in eq. (jlj) is over all triplets of integers j such that | j | < A and the limit A — > oo 
is implicit |l5|. This definition is equivalent to the analytic continuation of zeta-functions 
presented by Luscher . In eq. (JHJ) , L is the length of the spatial dimension in a cubically- 
symmetric lattice. The energy eigenvalue E n and its deviation from the sum of the rest 
masses of the particle, AE n , are related to the center-of-mass momentum p n , a solution of 
eq. ©, by 



AE n = E n - mi - m 2 = y pi + m\ + \j p\ + m\ - m 1 - m 2 
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where //12 is the reduced mass of the system. In the absence of interactions between the 
particles, |p cot <5| = 00, and the energy levels occur at momenta p = 27rj/L, corresponding 
to single-particle modes in a cubic cavity. Expanding eq. (j3J) about zero momenta, p ~ 0, 
one obtains the familiar relation 2 



AEn 



2ira 



1 + ci- + c 2 - 



(6) 



with 



1 



|j]<A 



Cl 



1 E t4 - 4A 



-2.837297 



j^o Ul 

and a is the scattering length, defined by 



c 2 
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tan 
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^2 : 14 



TT 



6.375183 , (7) 



(8) 



For the / = | ttK scattering length, 03/2, that we compute in this work, the difference 
between the exact solution to eq. (j3J) and the approximate solution in eq. (0) is much less 
than 1%. 



III. DETAILS OF THE LATTICE CALCULATION 

Our computation uses the mixed-action lattice QCD scheme developed by LHPC [16. 17] 
which places domain-wall valence quarks from a smeared-source on Nf = 2 + 1 asqtad- 
improved (lR Il9l | MILC configurations generated with rooted 3 staggered sea quarks [2^ 
that are hypercubic-smeared (HYP-smeared) |2£|, |29j, |3iJ, |3J| . In the generation of the MILC 
configurations, the strange-quark mass was fixed near its physical value, bm s = 0.050, (where 
b = 0.125 fm is the lattice spacing) determined by the mass of hadrons containing strange 
quarks. The two ligh t q uarks in the configurations are degenerate (isospin-symmetric) . As 
was shown by LHPC 0,01 j HYP-smearing allows for a significant reduction in the residual 
chiral symmetry breaking at a moderate extent L s = 16 of the extra dimension and domain- 
wall height M5 = 1.7. Using Dirichlet boundary conditions we reduced the original time 
extent of 64 down to 32. This allowed us to recycle propagators computed for the nucleon 
structure function calculations performed by LHPC. For bare domain-wall fermion masses 
we used the tuned values that match the staggered Goldstone pion to few-percent precision. 
For details of the matching see Refs. 0, Il7l |. The parameters used in the propagator 
calculation are summarized in Table |1J All propagator calculations were performed using 
the Chroma software suite 3^, 13 ■ 

As it is the difference in the energy between interacting mesons and non-interacting 
mesons that provides the scattering amplitude, we computed the one-pion correlation func- 
tion CV+(t), the one-kaon correlation function Cx+{t), and the kaon-pion correlation function 



2 We have chosen to use the "particle physics" definition of the scattering length, as opposed to the "nuclear 
physics" definition, which is opposite in sign. 



3 For recent discussions of the "legality" of the mixed-action and rooting procedures, see Ref. (20i I21L |2S 
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TABLE I: The parameters of the MILC gauge configurations and domain-wall propagators used in 
this work. The subscript I denotes light quark (up/down) where s denotes the strange quark. The 
superscript dwf denotes the bare quark mass for the domain wall fermion propagator calculation. 
The last column is number of configurations times number of sources per configuration. 



Ensemble 


bmi 


bm s 


i dw f 

bm l 


bmf»f 


10 3 x bm res a 


# of propagators 


2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


1.604 ± 0.038 


468x4 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


1.552 ±0.027 


658x4 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


1.239 ± 0.028 


486x3 


2064f21b681m030m050 


0.030 


0.050 


0.0478 


0.081 


0.982 ± 0.030 


564x6 



"Computed by the LHP collaboration. 



C n +K+{p, t), where t denotes the number of time-slices between the hadronic-sink and the 
hadronic-source, and p denotes the magnitude of the (equal and opposite) momentum of 
each meson. The single-pion correlation function is 

CV+W = $>•-(*, x)tt + (0,0)) , (9) 

X 

where the summation over x corresponds to summing over all the spatial lattice sites, thereby 
projecting onto the momentum p = state. The single-kaon correlation function has a 
similar form. The ir + K + correlation function that projects onto the s-wave state in the 
continuum limit is 

cw+GM) = £ E e * p ' &t " y) ^"(*.x)ir-(t l y) ^(0,0)^(0,0)) , (io) 
\p\=p x >y 

where, in eqs. © and (jIDj). 7r + (t,x) = u(t, x)j5(l(t, x) is an interpolating field for the 7r + , 
and K + (t,x) = u(t, x)7 5 s(t, x) is an interpolating field for the K + . In the relatively large 
lattice volumes that we are using, the energy difference between the interacting and non- 
interacting two-meson states is a small fraction of the total energy, which is dominated by 
the masses of the mesons. In order to extract this energy difference we formed the ratio of 
correlation functions, G w + K +(p, t), where 

<W(p,f) - - £ A « , (11) 



C 7T +(t)C K +(t) 



n=0 



and the arrow becomes an equality in the limit of an infinite number of gauge configurations. 
In G n +K+(p,t), some of the fluctuations that contribute to both the one- and two-meson 
correlation functions cancel, thereby improving the quality of the extraction of the energy 
difference beyond what we are able to achieve from an analysis of the individual correlation 
functions. 



IV. ANALYSIS AND CHIRAL EXTRAPOLATION 



A convenient way to present the data is with "effective scattering length" plots, simple 
variants of effective mass plots. The effective energy splitting is formed from the ratio of 
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correlation functions 

which in the limit of an infinite number of gauge configurations would become a constant at 
large times that is equal to the lowest energy of the interacting kaon and pion in the volume. 
At each time-slice, AE w + K +(t) is inserted into eqs.© and ©, or into eq. (JHJ), to give a 
scattering length at each time slice, a n + K +(t). It turns out to be more useful to consider 
the dimensionless quantity of the reduced mass times the scattering length, fi nK a n +x+, in 
our analysis, where the "effective reduced mass" is constructed from the effective 

mass of the single particle correlators. For each of the MILC ensembles that we analyze, the 
effective scattering lengths are shown in fig. [T] We found that the plateaus in [m w k i-k+k+ 
are better defined for smeared sinks than those from point sinks. 

The energy-shifts of the ir + K + groundstate can be extracted directly from the effective 
scattering length function, or from correlation functions, giving the same results. A single 
exponential function was fit by ^-minimization to the correlation functions, from which 
either the n + K + energy or single particle masses were determined. The central value and 
uncertainty of each parameter was determined by the jackknife procedure over the ensemble 
of configurations. The results of our lattice calculation of the decay constants, meson masses, 



TABLE II: Results from the lattice calculation. All errors are computed from jackknife. The 
uncertainty associated with the mOlO ensemble ttK energy shift and related quantities is dominated 
by the systematic error. The fitting ranges are shown in the square brackets. 



Ensemble 




m.K/U 




SE^K (MeV) 




T x 10 3 


m007 


2.000(17) 


3.980(25) 


1.332(10) 


11.89(81) [8-15] 


-0.1263(75) 


-10.1(9) 


mOlO 


2.337(11) 


3.958(16) 


1.469(07) 


11.40(50) [7-15] 


-0.155(40) 


-8(3) 


m020 


3.059(12) 


3.988(15) 


1.731(07) 


10.15(69) [10-15] 


-0.213(12) 


-5.59(34) 


m030 


3.484(10) 


4.004(12) 


1.869(05) 


10.06(54) [11-16] 


-0.267(12) 


-4.29(15) 



AE v+K+ (t) = log 



n + K + energy shifts and scattering lengths are tabulated in Table ITT1 The scattering lengths 
as a function of reduced mass are shown in fig. |2] 

As can be seen from fig. Q there is a large systematic error associated with the mOlO 
ensemble. There would appear to be two distinct plateaus. Rather than fitting to one of 
the plateaus, we chose to fit over a large range (7-15) which includes both plateaus and then 
assigned a systematic error which encompasses minima and maxima over the fit range as 
indicated by the effective scattering length plot. More statistics will have to be acquired on 
this ensemble before any conclusions can be drawn about this correlator. 

In SU(3) chiral perturbation theory 34[ 35|, 36|| at NLO, the expansion of the crossing even 
(a + ) and crossing odd (a~) scattering length times the reduced mass is known to be p, [D, 
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FIG. 1: The effective n + K + scattering length times the reduced mass, fi n jc a n +K+(t) as a function 
of time slice arising from smeared sinks. The solid black lines and shaded regions are fits with 1-a 
errors tabulated in Table The dashed lines on the mOlO ensemble plot are an estimate of a 
systematic error due to fitting. 
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FIG. 2: /jL^k 0-W+K+ vs - V"kk/ fir- The data points are the results of this lattice calculation, while 
the curve is the theoretical prediction at tree level in chiral perturbation theory The dark 

error bar is statistical, while the lighter error bar corresponds to the systematic error. The vertical 
dashed line denotes the physical pion and kaon masses. 
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where the counterterm L 7r ^-(A) is a renormalization scale, A, dependent linear combination 
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of the Gasser-Leutwyler counterterms 



L vK = 2Li + 2L 2 + L 3 - 2L 4 - y + 2L 6 + L 8 . (15) 

It is important to note that the expressions in eqs. (fTHj) and (fPfj) are written in terms of 
the full f n , and not the chiral limit value. The functions -yA NLO ' + ) (m^/A, mjf/A, m v /X) 
and „-/A, m^/A, m v /X) clearly depend upon the renormalization scale A. In the 

analysis that follows, it was found to be convenient to normalize the meson masses to /„., 
and therefore we can choose the renormalization scale to be A = f^ hys , and use the values of 
n^n/fn and mx/fn in Table ITT1 directly Deviations between the A = f n calculated on each 
lattice and A = /P hys are higher order in the chiral expansion. 

The I = \ and I — | scattering lengths are related to those in eqs. and (f*HJ) by 

a i/2 — a+ + 2a~ 

a 3/2 = a + - cT = a 7T+K + . (16) 

It is convenient to define the function F via a subtraction of the tree-level and one-loop 
contributions in order to isolate the counterterms, 

r ?£) s - JL [„„ K + 1 + X ^ - 2^ X ^°A ,(17) 

where we use the Gell-Mann-Okubo mass-relation among the mesons to determine the 77- 
mass, which we do not measure in this lattice calculation. At NLO this becomes 



r = L 5 (/f ys ) - 2 ^ L nK (ft yS ) ■ (18) 



It is clear that the dependence of F on and mx determines L5 and L n x and, in turn, 
allows an extraction of 03/2 and a\/2- The numerical values of F and their jackknife errors 
calculated on each ensemble of lattices are given in Table HH and are plotted in fig. El By 
fitting a straight line to the values of F as a function of mk/m^ the counterterms L 5 and 
L n x (renormalized at f% hjs ) can be determined. 

Ideally, one would fit to lattice data at the lightest accessible values of the quark masses 
in order to ensure convergence of the chiral expansion. While we only have four different 
quark masses in our data set, with pion masses ranging from m n ~ 290 MeV to 600 MeV, 
fitting all four data sets and then "pruning" the heaviest data set and refitting provides a 
useful measure of the convergence of the chiral expansion. Hence, in "fit A", we fit the data 
from all four lattice ensembles (m007, mOlO, m020 and m030), while in "fit B", we fit the 
data from the lightest three lattice ensembles (m007, mOlO and m020). 

With the limited data set presently at our disposal, it is not practical to fit to the NNLO 
expression for the scattering length. However, it is important to estimate the uncertainty 
in the values of the scattering lengths extrapolated to the physical point that is introduced 
by the truncation of the chiral expansion at NLO. In our work on fx/ fir we extracted 
a value of L 5 as it is the only NLO counterterm that contributes. The numerical value 
obtained is only perturbatively close to its true value, as it is contaminated by higher-order 
contributions. Therefore, by fixing the L 5 that appears in eq. (fTS|) to the value of L 5 extracted 
from fx/ fir, an estimate of the uncertainty in both L^k and in the extrapolated values 
of the scattering lengths due to the truncation of the chiral expansion can be estimated. 
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FIG. 3: r vs. mx/nT-Tr- The dark error bar on the data points is statistical, while the lighter error 
bar corresponds to the systematic error. The lines correspond to the four linear fits (A,B,C,D). 
The bars on the y axis represent the 1-a errors in the determinations of L5 = T(mx jfn-n = 0) as 
given in Table UTTI (At 95% confidence level, these determinations are in agreement.) 



TABLE III: Results of the NLO fits. The values of and m n ai/2 correspond to their ex- 

trapolated values at the physical point, where the error ellipses in the L^-L n K plane have been 
explored at 68% confidence level (see fig. @}. 

FIT L 5 x 10 3 L wK x 10 3 m n a 1/2 X 2 /dof 

A 3.83 ±0.49 3.55 ±0.20 -0.0607 ± 0.0025 0.1631 ± 0.0062 0.17 

B 2.94 ±0.07 3.27 ±0.02 -0.0620 ± 0.0004 0.1585 ± 0.0011 0.001 

C 5.65 ± 0.02±g;^| a 4.24 ±0.17 -0.0567 ± 0.0017 0.1731 ± 0.0017 0.84 

D 5.65 ± 0.02ig;£f a 4.16 ±0.18 -0.0574 ± 0.0016 0.1725 ± 0.0017 0.90 



Q Input from f K /U 37]. 



Specifically, we sampled L 5 from a Gaussian distribution for a range of fx/ fn values [37[ 
and then fit L^k using ^-minimization. We then generated a value of L n K from a normal 
distribution formed from its mean and standard error. This fit is denoted "fit C" , and the 
same fit but with the m030 data pruned is denoted "fit D" . The results of the four fits are 
given in Table 11111 and plotted in fig. El These fits lead to an extraction of 

L^k = 4.16 ±0.1818$ , (19) 
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2 I i I i I i I i I i I i I i 

1 2 3 4 5 6 7 

L 5 x 10 3 

FIG. 4: Error ellipses in the L^-L^k plane for the four fits (A,B,C,D) at 68% (dotted lines) and 
95% (solid lines) confidence level. 



and a prediction of the scattering lengths extrapolated to the physical point of 

m. a 3/2 = -0.0574 ±0.00161^058 

m v a l/2 = 0.1725 ± 0.0017±°;°i56 . (20) 

We have chosen to take the central values and statistical errors from fit D and have set the 
systematic error due to truncation of the chiral expansion by taking the range of the various 
quantities allowed by the four fits, including statistical and systematic errors. In fig. 0] we 
plot the 68% and 95% confidence- level error ellipses for the four fits given in Table ITTT1 in the 
L^L^k plane. In fig. El we plot the 95% confidence- level error ellipses associated with the 
four fits in the m w aM^-m^ a 3 / 2 plane 4 . For purposes of comparison we have included the 
current-algebra point |l0j on the plot as well as 1-a error ellipses from analyses based on 
fitting experimental data using x?T at NLO 0] and using Roy-Steiner equations 0. As 1-a 
error ellipses correspond to 39% confidence level, one should be careful in finding discrepancy 
between the various determinations of the scattering lengths. It would be interesting to see 
the NLO xPT and Roy-Steiner error ellipses at higher confidence levels. 



4 In Mathematica format, the 95% confidence- level error ellipses are: 
fit A: Ellipsoid[{0. 1631 ,-0.0607}, {0.0197,0.0007}, {{0.9283,0.3719}, {-0.3719,0.9283}}] 
fit B: Ellipsoid[{0. 1585,-0.0620}, {0.0076,0.0004}, {{0.9461,0.3239}, {-0.3239,0.9461}}] 
fit C: Ellipsoid[{0.1731,-0.0567}, {0.0042,0.0016}, {{0.7534,0.6576}, {-0.6576,0.7534}}] 
fit D: Ellipsoid[{0.1725,-0.0574},{0. 0046,0.0027}, {{0.7881, 0.6156}, {-0.6156,0. 7881}}] 
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FIG. 5: Error ellipses for the four fits (A,B,C,D) at 95% confidence level. (Note that these results 
are derived from lattice data on a single lattice spac ing of b = 0.125 fm.). The star corresponds to 
the current-algebra predictions (xPT p 2 ) from Ref. [l0(. We also display 1-er error ellipses from a 
XPT analysis at NLO (denoted xPT p A ) and from a fit using the Roy-Steiner equations [jj. 



Given how well our lattice data fit the NLO continuum %PT formulas, it would seem 
that the 0(b 2 ) discretization errors are comparable or smaller than the systematic error due 
to omitted 0{m? q ) effects in the chiral expansion. However, one should keep in mind that 
our determinations of, for instance, the low-energy constants L5 and L n K are subject to 
0(b 2 ) shifts. In contrast with the 7r + 7r + and K + K + scattering lengths, the mixed-action 



38, 39 . While 



quantity Amix makes an explicit contribution to the K + ir + scattering length 
this adds an additional unknown contribution to this process, a mixed-action xPT analysis 
of itK scattering, including lattice data from the fine MILC lattices (6 ~ 0.09 fm), will be 
able to address this source of systematic error quantitatively. We continue to search for the 
computational resources to accomplish this task. 



V. CONCLUSIONS 



In this paper we have computed the ir + K + scattering length in fully-dynamical lattice QCD 
at pion masses ranging between ~ 290 MeV and 600 MeV. We have used the continuum 
expressions for the scattering lengths in SU(3) chiral perturbation theory, together with 
lattice data for fx/ fir, to predict the physical I = 3/2 and / = 1/2 ttK scattering lengths 
with unprecedented accuracy. Naively one would expect that tt + K + scattering would give 
information about 1 = 3/2 scattering only. However, the lattice data, when combined 
with chiral perturbation theory, implies a constraint on / = 1/2 scattering as well. We 
anticipate that with improved statistics, together with calculations on lattices with smaller 
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lattice spacings, the theoretically-predicted regions for and can be further 

reduced beyond those shown in fig. |SJ These regions can then be compared with the expected 
measurements from K + tt~ atoms, to provide an exciting test of hadronic theory. 
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